Measurement-based quantum computer in the gapped ground state 

of a two-body Hamiltonian 
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We propose a scheme for a ground-code measurement-based quantum computer, which enjoys two 
major advantages. First, every logical qubit is encoded in the gapped degenerate ground subspace 
of a spin-1 chain with nearest-neighbor two-body interactions, so that it equips built-in robustness 
against noise. Second, computation is processed by single-spin measurements along multiple chains 
dynamically coupled on demand, so as to keep teleporting only logical information into a gap- 
protected ground state of the residual chains after the interactions with spins to be measured are 
turned off. We describe implementations using trapped atoms or polar molecules in an optical 
lattice, where the gap is expected to be as large as 0.2 kHz or 4.8 kHz respectively. 

PACS numbers: 03.67.Lx, 03.67.Pp, 42.50.Ex 
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Introduction. — Reliable quantum computers re- 
quire hardware with low error rates and sufficient re- 
sources to perform software-based error correction. One 
appealing approach to reduce the massive overhead for 
error correction is to process quantum information in 
the gapped ground states of some many-body interac- 
tion. This is the tactic used in topological quantum 
computation and adiabatic quantum computation. Yet, 
the hardware demands for the former are substantial, 
and the fault tolerance of the later, especially when re- 
stricted to two-body interactions, is unclear [1]. On the 
other hand, measurement-based quantum computation 
(MQC), in particular one-way computation on the 2D 
cluster state [2], runs by beginning with a highly entan- 
gled state dynamically generated from nearest-neighbor 
two-body interactions and performing computation by 
only single-qubit measurements and feed forward of their 
outcomes. However, its bare implementation may suffer 
decoherence of physical qubits waiting for their round of 
measurements in the far future, and that severely dam- 
ages a prominent capability to parallelize computation. 
Although its fault-tolerant method by error correction 
has been well established [3] , it is clearly advantageous if 
some gap-protection is provided on the hardware level. 

In this paper, we propose a ground-code measurement- 
based quantum computer (GMQC), which enjoys the two 
aforementioned advantages. GMQC is a conceptual ad- 
vance, since measurements generally create excitations 
in the system so that two desired properties, keeping the 
information in the ground state and processing the in- 
formation by measurements, are not seemingly compat- 
ible. We demand three properties of the ground state 
for GMQC. (i) There should be an energy gap to pe- 
nalize errors moving outside the computational ground 
subspace, and operators connecting logical states should 
be highly nonlocal. These should persist in the ther- 
modynamic limit to be scalable, (ii) The interactions 



should be preferably two-body. It is possible that in 
some ground subspace of H the effective interaction is 
if -local, but a demerit is the significantly reduced mag- 
nitude of the perturbative coupling, (iii) The interactions 
should be frustration-free, so that when every single spin 
is measured through computation, the remaining entan- 
gled spins to contain logical information can be set in the 
ground state of their Hamiltonian. 

We briefly refer to entangled resource states universal 
for MQC in the literature. First, the idea to use a ground 
state of the two-body Hamiltonian for universal MQC is 
seen in Ref. [4] , where the effective five-body Hamiltonian 
for the 2D cluster state is perturbativcly approximated 
in a low energy sector using ancillas. We do not resort 
to perturbation, so that not only is the gap much larger 
in practice, but also the resource state is exact in that 
we could approach unit fidelity as close as possible, by 
improving accuracy of analog simulation of our Hamilto- 
nian and its preparation to the ground state. Second, the 
novel use of tensor network states beyond the 2D cluster 
state [5, 6] has been quite motivating for us. But, our key 
idea of processing the logical information in degenerate 
ground states while maintaining the Hamiltonian on is 
incorporated for the first time in our paper. 

Scheme of a ground-code MQC. — In our GMQC, 
sketched in Fig. 1 , we adopt a hybrid approach where the 
logical two-qubit gates are implemented via dynamical 
couplings on demand as in the quantum circuit model, 
in addition to the standard MQC for the time evolution 
of each logical qubit. This is partly because there has not 
been known any exact gapped ground state of a two-body 
Hamiltonian which per se serves as a universal resource 
for the standard MQC. We utilize space-time resources 
in such a joint way that a "spatially" entangled resource 
is used to simulate the logical time evolution and "tem- 
poral" interaction is used to simulate the logical spatial 
interactions. Consequently, our GMQC exhibits some 
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FIG. 1: A schematic of a ground-code measurement-based 
quantum computer. Each logical qubit is encoded in every 
AKLT chain of N spin-1 particles with spin- 1/2 at the bound- 
aries and N + l, where spin particles interact with the near- 
est neighbors. The initial state of the chain is a unique ground 
state with a constant gap, and the measurement of the 0-th 
spin prepares the logical wire in \0 L } or \1 L ). Thereafter, com- 
putation is carried by measurements right to left, in process- 
ing the logical information in the two-fold degenerate, gapped 
ground state of each residual chain. Every single-qubit gate 
on the wire A is performed non-deterministically at the time 
step j by turning off the interaction PAj,A j+1 an d measuring 
the spin Aj in an appropriate basis. A non-deterministic two- 
qubit gate, CPHASEs.c, at the time step j is done by turning 
off both Pl jtBj + 1 and Pc j>Cj+1 , and then applying a dynam- 
ical coupling exp(i7r_ff m /x) followed by measurements of the 
spins Bj,Cj. Gate failures, which occur with bounded prob- 
abilities, are heralded, and act just as the logical identity by 
the adaption of subsequent measurement bases. The final 
readout is made by measuring the spin- 1/2 at TV + 1. 



new features in contrast with the conventional MQC. 

The basic Hamiltonian we consider is the ID Affleck- 
Kennedy-Lieb-Tasaki (AKLT) model [7], the chain 
of nearest-neighbor two-body interacting spins, H = 

J\Z?Jl] Ph + l + Po!l + Pn /2 N+iI With J > and Pfj+1 

the projector onto the spin-5 irreducible representation 
of the total spin for particles j and j + 1. Namely, 
P? j+1 = i(S, ■ S j+1 + .IS, • S^) 2 ) + V and pffl 



|(l6 + Sj ■ Sj>), where S,s are spin-1, 1/2 representa- 
tions of su(2). Without the boundary spins, the finite 
AKLT spin chain has a four-fold degeneracy correspond- 
ing to a total spin-0 (singlet) state and a triplet of spin-1 
states. The boundary terms P 3 / 2 project out the spin- 
1 components yielding a unique ground state such that 
H\Q) = 0. The gap AE persists in the thermodynamic 
limit, and has been estimated to be AE w 0.350J [8]. A 
key feature of \Q) is that it can serve, using single-spin 
measurements only, as a logical quantum wire which is 
capable of performing not only deterministic teleporta- 
tion as already remarked in Ref. [9] but also an arbitrary 
logical single-qubit operation similarly to one-way com- 
putation as shown below. 



uct state (MPS) representation [9, 10]. 
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(X, Y, Z are the 
Pauli matrices cr M (// = x,y,z)), and l^ - ) is the singlet 
(S = 0) located on the 0-th and + 1-th sites. This rep- 
resentation is helpful to see the action of local measure- 
ments to the "relative" state of unmeasured parts in sim- 
ulating unitary evolution of MQC [5, 6, 11]. Note that, in 
Refs. [5, 6], the ground state of a modified AKLT chain, 
different from the original, was considered to construct a 
resource state. However, their extension of the so-called 
byproduct operators into a non-Pauli finite group is less 
convenient, and it is unclear whether it still possesses a 
counterpart of the nonlocal string order utilized later. 

Universal computation by measurements. — 
Computation follows by measuring right to left (as in 
reading the Japanese comics) the single spins along each 
chain, where spins are indexed by increasing value mov- 
ing right to left, in accordance with the order of unitaries 
we simulate in the bra-ket notation. First we need to 
prepare the unique ground state \Q)® n of the n parallel 
decoupled ID AKLT chains. This can be done efficiently 
and deterministically along each chain, by either turning 
on H immersing the system in a reservoir and cooling it, 
or by making use of its MPS description to produce it 
via sequential unitaries (which scale linearly in N) [12] 
before turning on H. 

The initialization of every quantum logical wire is done 

3/2 

by first turning off the coupling P \ and measuring the 
rightmost 0-th spin-1/2 in the z basis. Because of the 



singlet configuration 



N+l/ 



the -1/2 and 1/2 out- 



come, denoted by |ln) and |0n) respectively, at the 0-th 
site 5-basis measurement induces the initialization of the 
wire (identified effectively with the preparation of the 
state at the N + 1-th site, since we will see that unitary 
actions accumulate on this degree of freedom) to \0 L ) 
and \l L ), respectively. With one boundary spin-1/2, the 
ground state is two-fold degenerate spanned by (lo\G) 
and (0n|{7), and computation takes place in this ground 
subspace after this initialization. 

Every logical single-qubit unitary operation is imple- 
mented by the single-spin measurement. The interaction 
Pj j+i is turned off before the local measurement of the 
j-th spin, to guarantee that the remaining system stays 
in the ground state. Since an arbitrary single-qubit op- 
eration is decomposed into three rotations around the 
logical Z and X axes with three Euler angles, wc show 
their measurement directions. The rotation R z (9) = 



3 



|0 L )(0 L | + e i9 \l L )(l L \ along the Z axis is applied by the 
single-spin measurement in an orthogonal basis, 

{|7f(f)» = {i((l±e- <tf )|l J -) + (lTe- <tf )|2 i )),|3 J )}. (2) 

If the outcome is either the first or the second (which 
turns out to occur always with probability 1 /3 for each) , 
we can apply R z (6) newly on the logical qubit with a 
byproduct operator X or XZ, respectively. If the out- 
come is the third, we interpret we have applied the "log- 
ical identity" with a byproduct Z. On the other hand, 



the rotation R x (6) = \+ L ){+ 1 
X axis, where |± L ) 



J6\ 



along the 



^(|0 L ) ± |1 L )), is applied by the 
single-spin measurement in another orthogonal basis, 

{bf (*)>} = {£((1 ± e ie )\2j) + (1 T e ie )\^)), II,-)}. (3) 

If the outcome is either the first or the second, we ap- 
ply R x (8) with a byproduct XZ or Z, respectively, and 
otherwise the logical identity with a byproduct X. 

Suppose we initialize the logical wire |0 L ) by the 0-th 
site —1/2 outcome (otherwise, we consider the wire is ini- 
tialized |0 L ) with the byproduct X from the beginning). 
According to Eq. (1), before the j-th spin is measured, 
we have a state = (1 | I\l>=j-i (7fc' \Q) S ivcn b Y 
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where the measurement directions of \"fk'{8k')) must be 
adapted from to — Qy when non-commuting byprod- 
ucts from previous measurements are propagated left 
through the current one, resulting in an accumulated 
byproduct operator T. 

We describe important properties of the 



Hamiltonian H(j) = J[Efe=/ P lk+i ^ J N,N+i. 
the measurement stage of the j-th spin. First. 
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residual 
through 
H(j) is 



gapped as before and is two-fold degenerate. Defining the 
string operators £ M (j) = e l7T s >= (£>cr^- +1 , we find that 
[^(j),H(j)] = whereas {£*(j), S 2 (j)}+ = 0. At each 
stage j, the pair Y, x ' z {j) forms a representation of su(2), 
and the degenerate ground states are only connected 
by nonlocal operators. For every single quantum wire, 
we utilize "time-dependent" logical encoding such that 
(mi^UMti)) = (0 N+1 \u^a» +i ru\0 N+1 ), where 
u is the total single-qubit rotation until the j — 1-th gate. 

Second, the logical state is not disturbed when turn- 
ing off the interaction coupling the bulk to the j-th spin. 
We can decouple with the time-dependent Hamiltonian 
H(j; t) = J{l-c{t))P'l J+l +H{j+l) , where c(i) is mono- 
tonically increasing in t E [0, 1] with c(0) = 0, c(l) = 1. 
Now Pjj+i does not commute with H(j + 1). However, 
the AKLT Hamiltonian has the property that the ground 
states also minimize its positive summands, i.e., it is 
frustration-free [7]. Hence, P gr (0) = P SI (t) Vf < 1 and 
g^p*l pg( t ) = Vi/,t, where P^ r (t) = \Q{t)){Q{t)\ are 



projectors onto the ground subspaces of H(j;t). Thus, 
turning off the end interaction term does not couple to 
excited states, and can be done in a constant time in- 
dependent of the system size. Even if there are some 
unwanted initial perturbations on iJ(j;i), the gap pro- 
vides robustness if performed adiabatically. The same 
argument applies to turning off the boundary terms. 

Third, notice that \ifj(j + 1)) of Eq. (4) can be written 
for the general outcome \tq) (r = 0, 1) of the 0-th site 
measurement as J2{a k }=i K+i) ■ ■ • l"A r )(?'o|[^ t ® 

rifetAr( a fe| M fc)]l*o,Af+i); where V is the form of 
■ R(8~k') due to the invariance of \^^) under the 
bilateral unitaries. But this is equivalent to the state 
obtained by beginning in the unique ground state of 



H(j + 1) + j+i> turn i n g °ff ^oj'+i' measuring the 0-th 
spin in the basis {V|t"o)> ^^l r o)}: with the result T^|ro) . 
This state is in the ground subspace of H(j + 1). 

Logical two-qubit operations are made dynamically by 
coupling two spin-1 particles in adjacent chains, say A 
and B, followed by their local measurements (equiva- 
lently, by a two-spin measurement). First P\. A , +i and 
^Bj B +i are turned off. We introduce the physical inter- 
action exp(ii/ lnt 7r/x) between spins Aj and Bj, where 
H** = X \S* A = 1}(S A . = 1| ® \S' Bj = l)(S* B . = 1|, 
and measure both spins Aj and Bj in the standard basis 
{|1), |2), |3)}. If the both outcomes are in either |1) or 
1 2), which occurs with overall probability (2/3) 2 = 4/9, 
we successfully apply the logical Controllcd-Phase gate 
CPHASE 4 ,b = U - 2|1^1^)(1^1^|. Notice that, in 
the span by \l Aj l Bj ), (1^2^^), \2 Aj l Bj ), and \2 Ai 2 Bj ) 7 
exp(iH int ir/ X ) acts as T = 1 4 - |(1 2 — X) ® (1 2 - X) 
and as the identity elsewhere. We see that this induces 
(a Aj \ ® (f3 Bj \r\M Aj ) ® |M Bj > = TCPHASEa^s with the 
byproduct T = XZ® XZ, XZ®X, X ® XZ, oxX®X 
in the aforementioned span. Otherwise (i.e., if at least 
one outcome is in |3)), we end up with applying the logi- 
cal identity with the byproduct (a Aj \M Aj ) ® (p Bj \M Bj ). 

To prove that computation is kept in the ground sub- 
space, imagine that we began in a separable state of two 
chains A and B initialized in |0 L ). After CPHASE is suc- 
cessfully applied at the stage j and byproducts for A and 
B are propagated, the joint state involves newly the op- 
erators (T A ®T b )\{1 4 + Z An+1 +Z Bn+1 -Z An+1 Z Bn+1 ), 
so that it is nothing but a superposition of logical states 
each of which is in the kernel of [H A (j + 1) + H B (j + 1)] 
and hence is in the ground subspace of the two chains. 

It can be verified that for any quantum circuit real- 
ized with our universal set of gates, the probability for 
each successful single-, two-qubit gate is constant at 2/3, 
4/9, respectively. We can efficiently perform the entire 
computation, by trying every gate until success, at the 
same time deterministically telcporting (by the standard 
basis measurement) other logical qubits to be spatially 
aligned for subsequent two-qubit gates. A remarkable 
new feature is that the general single-, two-qubit oper- 
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ations are probabilistic, while the logical identity (tele- 
portation) is essentially deterministic with the adaptive 
measurements. This variable computational depth orig- 
inates from the "correlated logical times" due to a two- 
point spatial correlation of the AKLT chain which decays 
as (— 1/3) I between sites j and f. 

At the end of computation, the joint state of left- 
most boundary spins is Y£/|0 L )® n , where T is the to- 
tal byproduct and U the target unitary operator. The 
final logical measurement, without loss of generality in 
the computational basis, is simulated by measuring in z 
these spins (after tclcporting the logical information if 
the chain is redundant). There, for example, the logi- 
cal |0 L ) outcome must be considered to occur if either 
s z N+1 = 1/2 and T docs not contain X, or s z v+1 = —1/2 
and T contains X. 

Physical implementations. — While the scheme of 
GMQC works independent of implementations, we sketch 
two physical realizations. In Refs.[13, 14], it was shown 
how to obtain the AKLT model (without the bound- 
ary terms) using tunneling induced collisions between 
bosonic atoms trapped one per site in a 3D optical lattice. 
The lattice beams are chosen such that the confinement 
is very strong along x, y directions but weaker along z, 
so that the system describes non-interacting ID chains. 
The spin states correspond to F = 1 hyperfine ground 
states of an alkali atom, and the effective coupling scales 
like J ~ t 2 /Us where t is the tunneling matrix element 
and Us (S = 0, 2) is a total spin dependent repulsive 
s-wave scattering energy in a common well. For realis- 
tic lattice parameters Us ~ 5 kHz and for t/Us = 0.3, 
we expect the gap AE m 0.2 kHz. A measurement on 
spin Aj is achievable by adiabatically expanding the lat- 
tice along z using the accordion technique [15] which de- 
creases t,Us but also allows the lattice wells within each 
wire to be better addressed. Then turning on optical 
tweezers [16] near the site Aj, one can shift the poten- 
tial depth at Aj so that tunneling to the neighboring 
wells is effectively zero. State dependent measurements 
can be done using spatially resolving microwave or Ra- 
man fields. To perform the CPHASE gate between Aj, Bj 
with wires separated along y, after turning off the inter- 
actions apply Raman pulses resonant at those locations 
that map \S Z = 1) to the first vibrational state of each 
well. Introducing a second lattice along y with twice the 
period of the first as in Ref. [14], the intensity can be 
adjusted so that the first excited vibrational states of 
nearest neighbor wells interact to generate a tunneling 
phase gate [17]. The second lattice intensity can be then 
tuned to turn off tunneling and the states mapped back 
to vibrational ground states to commence with measure- 
ment. To engineer the boundary terms, one could try to 
use different species on the boundaries with interactions 
satisfying P^, 2 = |(31 9 + 2sj • Sj, - Sf), where acts 
on in the basis {\S Z = 1), \ S Z = —1)} and zero elsewhere. 

An alternative is trapped spin-1 polar molecules with 



microwave-induced dipole-dipole interactions instead of 
tunneling along z. In Ref. [18], it was shown that the 
AKLT spin lattice model can be realized with J m 13.7 
kHz (AE f» 4.8 kHz). Boundaries can be loaded with 
spin- 1/2 species molecules using an additional optical lat- 
tice, and the necessary interactions could be designed us- 
ing spectroscopic resolvability for the microwave coupling 
fields between the different species. Interactions can be 
turned off by expanding the lattice along £ while tun- 
ing the microwave fields to keep the same interactions 
but with reduced strength. Then a rightmost molecule 
can be moved away from the bulk using optical tweez- 
ers and measured using state-dependent photoionization. 
The CPHASE gates can be generated between a pair of 
molecules in adjacent wires by outcoupling them using 
optical tweezers and using a microwave field to map 
\S Z = 1) -> \S Z = 0). Then H int is available using a sin- 
gle microwave field to induce the interaction (lg — S z2 )® 2 
[18], and measuring the spins in this new basis. In both 
implementations, the initial ground state can be prepared 
efficiently by adiabatically increasing the interaction of H 
from a configuration with antiparallcl spins induced by a 
staggered magnetic field as in Ref. [14]. 

Conclusion. — We have described a scheme to per- 
form MQC entirely within the gapped ground state 
of quantum many-body system with a two-body inter- 
action. For ground state protection, it is vital that 
the physical spin states are degenerate energy levels 
[19], meaning that local amplitude and phase errors are 
equally likely. The constant energy gap AE protects 
against noise whose spectral weight is smaller than it, 
providing a mechanism for anti-aging of the computa- 
tional resource. 
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